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Abstract The problem of the electromagnetic self-force can be studied in terms 
of a quadratic PT-symmetric Hamiltonian. Here, we apply a straightforward al¬ 
gebraic method to determine the regions of model-parameter space where the 
quantum-mechanical operator exhibits real spectrum. An alternative point of view 
consists of finding the values of the model parameters so that a symmetric operator 
supports bound states. 
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1 Introduction 

In a recent paper Bender and Gianfreda [1] discussed a model for the electromag¬ 
netic self-force on an oscillating charge particle. They showed that the problems 
reported by Englert [2] were due to the fact that the PT-symmetric Hamiltonian 
that describes the system exhibits broken PT symmetry. The addition of two in¬ 
teracting terms to that Hamiltonian enabled Bender and Gianfreda to analyse the 
regions of broken and unbroken PT symmetry. To this end the authors resorted 
to the approach of Rossignoli and Kowalski [3] in order to rewrite the quadratic 
Hamiltonian in diagonal form. The authors mentioned the possibility of using an 
alternative approach for the same purpose [3] which is known since long ago [5]. 
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The purpose of this paper is to stress the fact that the application of the latter 
algebraic method to quadratic Hamiltonians is extremely simple and straightfor¬ 
ward. In section [2] we derive the regular or adjoint matrix representation of the 
modified Hamiltonian operator for the electromagnetic self-force [1] and obtain the 
spectral frequencies in terms of the model parameters. These frequencies, which 
are roots of the characteristic polynomial for the matrix already mentioned, are 
suitable for determining the regions in parameter space where the PT symmetry is 
broken. Finally, in section [3] we summarize the main results and draw conclusions. 


2 The quantum-mechanical model 


From the pair of classical equations of motion proposed by Englert [2] , Bender and 
Gianfreda [1] derived the Hamiltonian function 


He 


PxPw — PyPz 2pzPw WPy -h zpx 

mr rriT^ 2 


mzw , 

—2 -f kxy, 


( 1 ) 


where q = x,y,z,w are suitable coordinates and pq the corresponding conjugate 
momenta. The quantum-mechanical version of this operator is PT-symmetric but 
its eigenvalues are not real because the PT symmetry is broken for all m, r, k. To 
overcome this difficulty Bender and Gianfreda added two terms and obtained the 
modified Hamiltonian operator 


B {wpz - ZPw) 2pzPw ^ PxPw - PyPz 
mr mr'^ mr 

A (x^ + 1/2) 

^ 2 
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~2 ^ 2 
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( 2 ) 


where the coordinates and momenta satisfy the quantum-mechanical commutation 
relations [q,Pq] = i- In this way those authors were able to show that the resulting 
PT-symmetric Hamiltonian exhibits real spectrum for suitable values of A and B. 
In particular, the PT symmetry is broken for the case A = B = 0 that leads to 
He. 

The simplest way of determining the regions of broken PT symmetry is to 
convert H into a suitable diagonal form. Bender and Gianfreda resorted to the 
approach proposed by Rosignoli and Kowalski [3] and in what follows we apply 
a well known algebraic approach[F96] proposed recently by Fernandez [3] for the 
treatment of a simple model for optical resonators [6]. 
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The algebraic approach is based on the construction of the adjoint or regular 
matrix representation of H in the operator basis {Oi, O 2 , ■ ■ •, Os} = y, z, w,px,Py,pz,Pw}- 
The matrix elements H^j are given by the coefficients of the commutator relations 


i=i 


A straightforward calculation leads to 

/n n n 0 
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(4) 


The spectral frequencies of the problem are the eigenvalues of this matrix. The 


secular equation |H — AI| = 0 (where I is the 8x8 identity matrix) yields the 
characteristic polynomial 


/ 2 2, o 2 I 2\ / 2 2*3 , ^2 / 2 
(m T ^ — B + m ) [m r ^ ^ [m 

where ^ One of the roots is 


+ C (2AB - 2km) - A^ + k^) = 0, 

(5) 


C = 




and the remaining three ones are solutions to the cubic equation 


( 6 ) 


^2^2^3 ^^2 _ ^2 2 (AB -km)(-A^+k^= 0. (7) 

When the spectrum is real the frequencies are real and the four roots are positive. 
Thus the region in parameter space where the spectrum is real is determined by 
the set of equations A^ > k^, AB > fcm}. 

Bender and Gianfreda [T] obtained exactly the same equations m and 0 by 
means of an alternative approach based on the creation and annihilation operators 
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[3]- In our opinion the adjoint or regular matrix representation of the Hamiltonian 
operator provides a more straightforward route. In addition to it, the algebraic 
approach is somewhat more general because it applies to any Hamiltonian that 
satisfies the commutator relations m 0- 


3 Conclusions 

Any Hamiltonian that is a quadratic function of the coordinates and momenta 
or of the creation and annihilation operators is exactly solvable. One way of ob¬ 
taining its spectrum is to convert it into a diagonal form by means of a linear 
combination of the relevant dynamical variables [3] . The algebraic method [ills] is 
a particularly simple approach because the construction of the adjoint or regular 
matrix representation is straightforward. It only requires the calculation of simple 
commutation relations like ([3]). 

It is also worth mentioning that the quantum-mechanical versions of He and 
H are symmetric operators. They satisfy {f\H\g) = {g\H\f)* for any pair of 
square integrable functions /(q) and g{c^. If V' is a normalizable eigenfunction of 
H [Hip = Et/j, {t/j] Ip) < oo) then E is real. Therefore, in the present case, looking for 
the regions in parameter space of unbroken PT symmetry is equivalent to looking 
for the regions where H supports bound states (as Bender and Gianfreda [1] did 
explicitly for the ground state). The same situation emerged in the case of the 
optical resonators mentioned above HE]. 
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